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We present evidence for Rydberg-Rydberg interaction in a gas of rubidium atoms above room 
temperature. Rabi oscillations on the nanosecond timescale to different Rydberg states are investi- 
gated in a vapor cell experiment. Analyzing the atomic time evolution and comparing to a dephasing 
model we find a scaling with the Rydberg quantum number n that is consistent with van der Waals 
interaction. Our results show that the interaction strength can be larger than the kinetic energy 
scale (Doppler width) which is the requirement for realization of thermal quantum devices in the 
GHz regime. 
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Loosely bound electrons in highlying Rydberg states 
are giving rise to long range dipolar and van der Waals 
interactions [TJ [2] . The coherent control of this interac- 
tion allows for engineering of quantum correlated states 
[31 2] . In ultracold systems where an interaction strength 
on the order of several MHz is sufficient, the van der 
Waals interaction has lead to the observation of Rydberg 
blockade E] and dephasing [5] . In such systems the 
interaction has been exploited for quantum logical oper- 
ations [HI QH] and the creation of non-classical light states 
[TT] . So far the observation of this interaction has been 
limited to the ultracold domain. In very early related ex- 
periments density-dependent line broadening effects on 
Rydberg lines have been studied (T2j and interaction ef- 
fects in thermal vapor involving excited but non Rydberg 
states have been investigated [T3J El] • 

We present evidence for van der Waals-type inter- 
atomic interaction energies in the GHz regime between 
Rydberg-excited alkali atoms in thermal vapor. Using a 
pulsed laser excitation, we are able to drive Rabi oscil- 
lations on the nanosecond timescale to a Rydberg state 
[TBI and are therefore faster than the coherence time lim- 
itation given by the Doppler width. 

We investigate the dephasing of these oscillations for 
different atomic densities and Rydberg states. For a fixed 
Rydberg state, we see a linear growth of the dephasing 
rate with density. Through a systematic study of vari- 
ous Rydberg states we have found that the scaling of this 
growth with the principal quantum number n is consis- 
tent with van der Waals-interaction. 

We excite 85 Rb atoms to a Rydberg S-state with an 
off- resonant two-photon excitation (Fig. [TJi). The upper 
transition is addressed by a pulsed laser in order to pro- 
vide sufficient intensity for driving GHz Rabi oscillations 
despite the small transition dipole matrix element. The 
peak Rabi frequencies of the pulse on the 5P — nS transi- 
tion are in the range of f2biue/27r ~ 3.3 GHz to 3.7 GHz. 
The 780 nm laser addressing the ground state (fi re d/27r ~ 
750 MHz) is blue detuned by A/2tt ~ 2.3 GHz with re- 
spect to the intermediate state |15| . In this configuration 



a major part of the population oscillates directly between 
the ground and Rydberg state. In particular this means 
that changes in the absorption/emission of the 780 nm 
light are mainly caused by these population transfers and 
not due to excitation of the intermediate state. 

A crude estimate of the associated effective two-photon 
Rabi frequency is given by £l e ff/2ir w Oblue^red/(2A • 
27r) = 550 MHz. Therefore, the coherent dynamics is 
faster than the dephasing due to atomic motion expressed 
by the Doppler width of the two-photon transition of 
< 370MHz for our experimental parameters. Neverthe- 
less it has to be noted that the intermediate level cannot 
be completely neglected as is done usually in cold atom 
experiments, as the Rabi frequencies are comparable to 
the detuning. We will therefore model the experiment 
using three electronic states. The exact detuning for 
the 480 nm pulse is set to the point of resonant (slow- 
est) Rabi oscillation, which is determined by a detuning 
scan (Fig. [TJ;) . Note that the two-photon resonance expe- 
riences a time-dependent Autler-Townes shift to smaller 
frequencies due to the strong laser coupling. 

The experimental setup is essentially the same as dis- 
cussed in [IB]. The two laser beams are overlapped in 
a 5 mm vapor cell in almost counter-propagation con- 
figuration (Fig. [lb). The Rabi oscillations during the 
pulse translate to oscillations in the transmitted inten- 
sity of the 780 nm laser which is subsequently recorded 
by a fast photodetector. In order to prevent a washing- 
out of the oscillations due to a distribution of different 
Rabi frequencies, we use imaging systems before and af- 
ter the cell to ensure detection from a homogeneously 
illuminated region. 

We have been investigating the density dependence 
of the Rabi oscillations for densities N g in the range of 
10 11 — 10 13 cm" 3 and for different Rydberg S-states. As 
an example the data for n = 37 is shown in Fig. [2] Note 
that the curves have been normalized as the original sig- 
nal amplitude is approximately proportional to the den- 
sity. We observe an increasing damping in the oscillations 
with higher density which points towards a dephasing or 
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Fig. 1. Schema of the experiment, (a) Relevant level scheme 
and excitation wavelengths. We use a two-photon-excitation 
scheme for addressing a Rydberg S-state via the D2-line. The 
excitation is performed off-resonantly with a (blue) detuning 
A/27T of ~ 2.3 GHz with respect to the intermediate state, (b) 
Optical setup. The two beams are overlapped in the vapor 
cell in almost counter-propagation configuration. We are de- 
tecting the transmission of the 780 nm laser in which the Rabi 
oscillations manifest themselves during the excitation pulse, 
(c) Oscillation signals of a detuning scan of the 480 nm laser. 
We fix the frequency position for the analysis to the point of 
slowest oscillation (dashed line). 



loss mechanism of the Rydberg state. 

Ionization and Rydberg-ground state collisions do not 
give a significant contribution on the timescale of our 
experiment [IT]. Also incoherent plasma formation in 
Rydberg gases has been reported to be on a time scale 
of ~ 100 ns or larger (THJ Qjj] and can be neglected. 
In order to exclude broadening effects on single atoms 
like e.g. perturbation by electric or electromagnetic 
fields or collisions with a background gas, we checked for 
narrow EIT (electromagnetically induced transparency) 
lines in a CW Rydberg-EIT experiment 20 with low- 
power lasers. The observed line widths were less than 
20 MHz. As the Rydberg population in such an EIT ex- 
periment is very small (here < 0.1%), Rydberg-Rydberg 
interactions are not expected to contribute significantly 
to the broadening of the EIT line [211 E2]- Having ex- 
cluded single atom effects, we conclude that the dephas- 
ing can only be explained by an interatomic effect that 
is proportional to the Rydberg population. 

In order to quantify the damping behavior we employ 
a single atom 3-level model for the coherent dynamics 
[T5] , Besides the measured pulse shape as an input, the 
model includes two additional terms that describe a de- 
phasing of the Rydberg state. A constant dephasing term 
Toffset accounts for an offset dephasing that occurs in our 
data. This dephasing varies with the alignment of our 
pulse amplifier setup which has to be changed for differ- 
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Fig. 2. Rabi oscillations for different atomic densities. Os- 
cillations in the transmission of the 780 nm laser for different 
atomic densities are shown here for the Rydberg state 37S. 
The numbers on the right-hand side give the atomic ground 
state density in 10 12 cm~ 3 . The experimental data (black 
line) is averaged over 20 shots and normalized in amplitude. 
A curve derived from a 3-level-atom model (red line) is fitted 
to the data. The only fit parameters are the dephasing rate of 
the Rydberg state and the signal amplitude. The gray-shaded 
curve in the background depicts the measured temporal shape 
of the excitation pulse intensity. 



ent principle quantum numbers. However for a series of 
measurements at a particular principle quantum number 
it is constant. We attribute it - besides atomic motion - to 
phase noise and imperfections in the spatial homogeneity 
of the pulse. For the model r o ff sct is chosen such that it 
accounts for all remaining dephasing when extrapolating 
to zero atomic density. It is on the order of few hundred 
MHz. The density-dependent part of the dephasing is 
described by a term p 33 • r p33 that is proportional to the 
Rydberg population P33. 

The dephasing rates r p33 are extracted from the exper- 
imental data by fitting curves obtained from this model 
(Fig. [2j red lines) . The only fit parameters are r p33 and 
the signal amplitude. For each Rydberg state we find a 
linear increase in the dephasing rate with atomic density 
Ag (Fig. [3^i). The corresponding slope is strongly increas- 
ing with the principal quantum number of the Rydberg 
state. Note that the dephasing model only fits well to 
experimental signals below a certain density (depending 
on the Rydberg state). For higher densities the model 
does not accurately describe the data anymore such that 
it is not possible to extract values for dephasing rates. 
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Fig. 4. Slope of the dephasing vs. principal quantum number. 
Double- logarithmic plot of the dephasing slopes dT P33 /dN g 
vs. the effective principal quantum number n*. The fit of 
a power law oc n, p yields an exponent of p = 5.58 ± 0.48 
(solid line). This is in excellent agreement with an rc* 11 ' 2 
van der Waals scaling. Fitting an n* 4 scaling curve to the 
data (dashed line), it can be seen that resonant dipole-dipole 
interaction is not consistent with the data. The inset shows 
the dephasing at the van der Waals critical density. For N g = 
iVcrit all examined Rydberg states show the same dephasing. 



Fig. 3. Dephasing parameters vs. density, (a) The dephas- 
ing parameter of the Rydberg state r P33 corresponding to 
the Rydberg-Rydberg interaction is shown as a function of 
the atomic ground state density N s . For each Rydberg state 
the data is fitted by a linear curve, (b) When scaling the 
density-axis with the critical density N CT it all data points col- 
lapse onto a single line. -/V cr it is determined by the van der 
Waals-blockade radius for each Rydberg state individually. 
The dashed line is a linear fit to all data points. 

Whether these effects are due to collective many-body 
behavior is subject to further study. The 22S-state shows 
no increase in the dephasing rate within the margins of 
error which confirms that our offset dephasing is indeed 
constant during one series of measurements. 

In order to obtain a unified behavior we define a criti- 
cal density N crit that corresponds to the interatomic dis- 
tance at which the van der Waals-shift matches the ex- 
citation bandwidth O e g, which is kept constant for all 
experiments: 

htleS = — ~~c = Cq ■ N CTit 2 . 

Rescaling the atomic density in units of iV cr i t for each 
principal quantum number, we find that all dephasing 
values collapse onto a single line (Fig. [3]d). 

The van der Waals-shift is expected to behave like 

where = n — 6$ is the effective principal quantum 
number accounting for the quantum defect 5s ~ 3.13 |23j . 
Accordingly the critical density scales as N CT it oc rt^T 11 ! 2 . 



We therefore would expect the slopes of the dephasing 
rates to behave like dT P33 / dN g oc n* 11 / 2 . Fitting a power 
law oc nj 1 to the slopes for the different effective principal 
quantum numbers (Fig. H), we obtain an exponent of 
p = 5.58±0.48 in excellent agreement with van der Waals 
scaling. 

At the interatomic distances of ~ 1 /im present in the 
experiment the interaction potentials of different pair 
states can potentially be within the excitation bandwidth 
|24) . Hence, the perturbative limit which is van der Waals 
interaction might not necessarily be valid. If a second 
pair state (e.g. n'Pn"P) with an allowed dipole tran- 
sition to the respective nSnS state is within the excita- 
tion bandwidth, the interaction energy will also include a 
C^/r 3 term and thus exhibit n* -scaling [25]. This scal- 
ing, however, clearly lies outside the margins of error of 
the fit and can thus be ruled out. Therefore, we conclude 
that for the nSnS states examined the van der Waals de- 
scription remains valid for interparticle distances down 
to at least 1 /mi. 

Looking again at the dephasing rates at the van der 
Waals critical density, we can see that all Rydberg states 
dephase with the same rate of ~ 800 MHz • p 33 (Fig. |4j 
inset). Note that the definition of iV cr jt contains a some- 
what arbitrary constant factor which has, however, no 
influence on the n-scaling. Therefore, the value at -/V cr j t 
is reflecting only the order of magnitude of the interaction 
energy at the blockaded density. 

At higher densities (N g w 2N CV i t ) we have measured 
dephasing rates of up to ~ 2 GHz (for fully excited 
Rydberg atoms). This value clearly lies outside of the 
Doppler width of the sample. This shows that the un- 
derlying van der Waals interaction strength can be large 
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enough to support a blockade effect in a thermal gas. 

In summary we have found evidence of Rydberg- 
Rydberg interaction in the coherent evolution of atoms 
in thermal vapor. With a systematic study of the scaling 
for different Rydberg quantum numbers we could identify 
the interaction to be of van der Waals type. We have seen 
that the strength of the observed van der Waals interac- 
tion can be larger than the Doppler width. This is the 
basic requirement for the observations of a blockade effect 
and coherent collective behavior in a thermal ensemble. 
Hence these results suggest room temperature alkali va- 
pors as promising candidates for quantum devices. For 
the direct observation of a blockade effect further studies 
in a small spatial volume on the order of the interaction 
range [26] will be necessary. 

We acknowledge insightful discussions with J. Shaf- 
fer on the details of the molecular potentials and thank 
K. Rzazewski and H. Kiibler for helpful suggestions. The 
work is supported by the ERC under contract number 
267100, BMBF within QuOReP (Project 01BQ1013), the 
EU project MALICIA, and contract research 'Interna- 
tionale Spitzenforschung IF of the Baden- Wurttcmberg 
Stiftung. T.B. acknowledges support from the Carl Zeiss 
Foundation. B.H. acknowledges support from Studien- 
stiftung des deutschen Volkes. 



[13 



* These authors contributed equally to this work. 
' t.pfau@physik.uni-stuttgart.de 
[1] T. F. Gallagher, Reports on Progress in Physics 51, 143 
| (1988)| " 
[2] M. Saffman, T. G. Walker, and K. M0lmer, |Rev. Mod. 



Phys. 82, 231 3 (2010) 
[3] D. Jaksch, J. I. Cirac, P. Zoller, S. L. Rolston, R. Cote, 

and M. D. Lukin, |Phys. Rev. Lett. 85, 2208 (2000)| 
[4] M. D. Lukin, M. Fleischhauer, R. C ote, L. M~Duan, 
D. Jaksch, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 
87, 037901 (2001)| 
[5] D. Tong, S. M. Farooqi, J. Stanojevic, S. Kris hnan, Y. P. 
Zhang, R. Cote, E. E. Ey ler, and P. L. Gould, |Phys. Rev.| 
Lett. 93, 063001 (2004)] 
[6] R. Heidemann, U. Raitzsch, V. Bendkowsky, B. Butscher, 
R. Low, L. Sa ntos, and T. Pfau, |Phys. Rev. Lett. 99, | 
| 163601 (2007)| 
[7] K. Singer, M. Reetz-Lamour, T. Amthor, L. G. Marcassa, 
and M. Weidemuller, |Phys. Rev. Lett. 93, 163001 (2004)] 
[8] T. A. Johnson, E. Urban, T. Henage, L. Isenhower, D. D. 
Yavuz, T. G. Walke r, and M. Saffman, |Phys. Rev. Lett.| 
| 100, 113003 (2008)| 
[9] L. Isenhower, E. Urban, X. L. Zhang, A. T. Gill, 
T. Henage, T. A. Johnson, T. G. Walker, and 
M. Saffman, |Phys. Rev. Lett. 104, 010503 (2010)| 



[10] T. Wilk, A. Gaetan, C. Evellin, J. Wolters, Y. Miroshny- 
chenko, P. Grangier, and A. Browaeys, |Phys." Rev. Lett. 
I 104, 010502 (2010)| 

[11] Y. O. Dudin and A. Kuzmich, |Science 336, 887 (2012) 
J. M. Raimond, G. Vitrant, and S. Haroche, J. Phys. B 



12 



At. Mol. Phys. 14, L655 (1981)| 



14 



[15 
[16 



[is: 

[19 
[20 



[21 



[22 



F. Shen, J. Gao, A. A. Senin, C. J. Zhu, J. R. Allen, 
Z. H. Lu, Y. X iao, and J. G. Eden, |Phys. Rev. Lett. 99,| 
143201 (2007)| 

J. Keaveney, A. Sargsyan, U. Krohn, I. G. Hughes, 
D. Sarkisyan, and C. S. Adams, Phys. Rev. Lett. 108, 
173601 (20i2)| 

See supplemental material for further details on the ex- 
perimental setup, 3-level model and dephasing mecha- 
nism. 

B. Huber, T. Baluktsian, M. Schlagmiiller, A. Kolle, 
H. Kiibler, R. Low, and T. Pfau, |Phys. Rev. Lett. 107,| 
24300 H20TT)| 

D. C. Th ompson, E. Weinberger, G.-X . Xu, and B. P. 
Stoicheff, |Phys. Rev. A 35, 690 (1987)] 

G. Vitrant, J. M. Raimond, M. Gross, and S. Haroche, 
|J. Phys. B: At. Mol Phys. 15, L49 (1982) 



M. P. Robinson, B. L. Tolra, M. W. Noel, T. F. Gal- 
lagher, and P. Pillet, |Phys. Rev. Lett. 85, 4466 ( 2000) 
A. K. Mohapatra, T. R. Jackson, and C. S. Adams, Phys. 
Rev. Lett. 98, 113003 (2007)] 

H. Schempp, G. Giinter, C. S. Hofmann, C. Giese, S. D. 
Saliba, B. D. DePaola, T. Amthor, M . Weidemiiller, 
S. Sevi ncli, and T. Pohl, |Phys. Rev. Lett. 104,T7 3602 
(2010)| 

J. D. Pritchard, D. Maxwell, A. Gauguet, K. J. Weath- 
erill, M. P. A. Jones, and C. S. Adams, Phys. Rev. Lett. 
105, 19360 3 (20107] 

W. Li, I. Mourachko, M. W. Noe l, and T. F. Gallagher, 
|Phys. Rev. A 67, 052502 (2003)| 

J. S. Cabral, J. M. Kondo, L. F. Gonalves, V. A. 
Nascimento, L. G. Marcassa, D. Booth, J. Tallant, 
A. Schwettmann, K. R. Overstreet, J. Sedlacek, and J. P. 
Shaffer, |J. Phys. B: At. Mol. Phys. 44, 184007 (2011)| 
A. Reinhard, T! C~! Liebisch, R Knuffman, and 

G. Raithel, [Phys. Rev. A 75, 032712 (2"007)] 

H. Kiibler, J. P. Shaffer, T. Baluktsian, R. Low, and 
T. Pfau, |Nat. Photon. 4, 112 (2010)| 

A. Schwettmann, C. McGuffey, S. Chauhan, K. R. Over- 
street, and J. P. Shaffer, [Appl. Opt. 46, 1310 (200771 
M. Fleischhauer, A. Imamoglu, and J. P. Marangos, Rev. 



Mod. Phys. 77, 633 (2005) 



5 



SUPPLEMENTAL MATERIAL 



Experimental details. 

The ground state laser at 780 nm is blue detuned by 
~ 2.3 GHz with respect to the center of gravity of the 
5P 3 / 2 state. This is achieved by locking the laser onto 
a line of the neighboring 87 Rb isotope (5Si/ 2 ,-F = 1 — > 
5P 3 / 2 ) transition. The 480 nm laser pulse driving the 
upper transition is created with a seeded dye amplifier 
setup in 4-pass configuration similar to the one described 
in [27 . The transmitted light at 780 nm is detected with 
an AC-coupled avalanche photo diode with a lower cutoff 
at ~ 1 MHz. This means only fast changes in the signal 
are detected. The beam diameters at the position of the 
cell are 0hiuc = 185 pm and rc d = 220 pm. Further 
details about the optical setup can be found in previous 
work Hi . 

A reservoir containing a droplet of rubidium is at- 
tached to the 5 mm vapor cell and can be temperature- 
controlled. The atomic vapor density inside the cell can 
be adjusted by changing this temperature. The cell itself 
is placed in an oven and kept at a constant temperature 
of T ~ 130° C in order to prevent condensation of atoms 
on the glass surface. The exact atomic density in the 
cell is determined by absorption measurements. 



3-level model. 

As shown in previous work |16| , we model the time evo- 
lution of our system with a single atom three-level system 
with a density matrix approach using the Liouville-von 
Neumann equation |28j 



dp 
Of 



L (p) + L deph (p) 



H is the Hamiltonian containing the light couplings and 
L (pi) contains the decays of the levels due to natural 
linewidths. In order to account for the observed changes 
in the signal with increasing density we have introduced 
an additional dephasing of the Rydberg state of the form 



-^deph (p) — (r p33 P33 + r offset) 



-5031 



■§p32 




The dephasing term T P33 p S3 is proportional to the Ryd- 
berg population and thus describes Rydberg-Rydberg 
interaction. r o ff sct accounts for a constant offset dephas- 
ing that we attribute to experimental imperfections. 
The dephasings only act on the coherences of the density 
matrix such that no population transfer is induced 
by them. For the calculations we have neglected the 
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Fig. 5. Model for the dephasing behavior. Each atom experi- 
ences a level shift given by the van der Waals interaction po- 
tential (black curve). The distance distribution of the atoms 
is shown in the bottom part for a density of N = 0.25N cr i t 
(blue line) and N = 0.5N cr i t (red dashed line). According to 
the 1 /r 6 -potential this translates to a distribution of Rydberg 
level shifts shown on the upper left part (blue and dashed red 
curves, respectively). The excitation bandwidth (gray-shaded 
curve), which is mainly determined by the effective Rabi fre- 
quency fi, gives a feeling for which level shifts contribute to 
the signal. The critical density is defined such that the van der 
Waals shift at the corresponding mean distance r cr it matches 
the excitation bandwidth (black dashed lines) . 



velocity distribution of the atoms as the dynamics is 
dominated by the laser detunings (Aooppicr *C ALascr)- 
We have verified that the results with and without 
Doppler distribution differ by less than ~ 5%. The 
time dependence of the blue Rabi frequency f^biuc wa s 
extracted from the real temporal envelope of the laser 
pulse. A digital frequency filter was applied to the 
results to mimic the finite bandwidth (~ 1GHz) of our 
detector. 



Dephasing mechanism. 

In order to understand the dephasing mechanism more 
quantitatively we developed the following simple model. 
The model is based on an ensemble of atoms, in which 
the time-evolution of each atom is described by single- 
atom dynamics only. The Rydberg state of each atom 



experiences a van der Waals shift given by the distance 
to its nearest neighbor (Fig. [5]) . The distribution of inter- 
atomic distances in a thermal vapor is then mapped to a 
distribution of interaction strengths. Assuming a nearest 
neighbor distribution of w(r) = 4Trr 2 Nexp (— ^r 3 iV) for 
a certain atomic density AT, we calculate the Rabi oscilla- 
tions in the signal for each level shift and average then ac- 
cording to the distribution of interaction strengths. The 
resulting signal is evaluated in the same way as the exper- 
imental data by fitting curves from the dephasing model. 
The model reproduces the linear density dependence of 
the dephasing rates for densities N up to ~ 0.5iVcrit 
(Fig. [5j inset) . The assumption in the model is that each 
atom experiences the full shift of the interaction which 
implies that each nearest neighbor is fully excited to the 
Rydberg state. In the real experiment, however, only an 
average Rydberg population of ,033 ~ 20% — 25% (aver- 
aged over time) is achieved. This is consistent with the 
experimental data points ranging until N s ~ 2iV cr i t . Also 
the magnitude of the dephasing rate is in good agreement 
with the experimental results. 

Furthermore it becomes clear that there is a direct con- 
nection between the excitation bandwidth and the de- 
phasing rate as only atoms whose shifts lie within the 
bandwidth can contribute to the signal. The single atom 
dynamics in the model is justified as long as the frac- 
tion of blockaded atoms is small enough as they do not 
contribute to the signal then. For a significant number 
of blockaded atoms, however, collective dynamics is ex- 
pected to appear [6] and the model breaks down. 



